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Abstract 

We test whether or not the r lepton manifests the same couphngs as the /U lepton by investigating 
the relative decay rates in purely leptonic D'^ meson decays. Specifically, we place the first upper 
limit on the ratio R = T {D~^ -^ t~^v) /T {D~^ -^ f-^^)- We use 281 pb^^ of data accumulated at 
the V'(3770) resonance with the CLEO-c detector, to determine B{D+ -^ t^u) < 2.1 x 10"^ at 90% 
confidence level (C. L.). The ratio of R to the Standard Model expectation of 2.65 then is <1.8 at 

% C. L., consistent with the prediction of lepton universality. 



I. INTRODUCTION 

The Standard Model decay diagram for D^ -^ i'^u is shown in Fig. 1. The decay rate is 
given by [1] 




where Md+ is the D~^ mass, mi is the mass of the charged final state lepton, Vcd is a Cabibbo- 
Kobayashi-Maskawa matrix element with a value we take equal to 0.225 [2], and Gp is the 
Fermi coupling constant. 
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FIG. 1: The decay diagram for D^ ^ i'^iy. 

The decay is helicity suppressed because the virtual W~^ is a spin-1 particle, and the 
final state consists of a naturally left-handed spin- 1/2 neutrino and a naturally right-handed 
spin-1 /2 anti-lepton that have equal energies and opposite momenta. The ratio of decay 
rates for any two different leptons is then fixed by well-known masses. For example, for t^u 
to fJ,^!^, the expected ratio is 

,2 



mt+ 1 



TT 



Any deviation from this formula would be a manifestation of physics beyond the Standard 
Model. This could occur if any other charged intermediate boson existed that coupled to 
leptons differently than mass-squared. Then the couplings would be different for muons and 
r's. This would be a manifest violation of lepton universality, which has identical couplings 
of the muon, the tau, and the electron to the gauge bosons (7, Z^ and W^) [3]. (We note 
that in some models of supersymmetry the charged Higgs boson couples as mass-squared to 
the leptons and therefore its presence would not cause a deviation from Eq. 2 [4].) Using 
measured masses [5], this expression yields a value of 2.65 with a negligibly small error. 

We have already reported [6] B{D+ -^ fi+u) = (4.40 ± OMtoHl) x 10"^, and established 
an upper limit of B{D^ -^ e^v) < 2.4 x 10^^. It remains to measure or limit t^v, which is 
the subject of this paper. We note, for reference, that the predicted relative widths in the 
Standard Model are 2.65 : 1 : 2.3 x 10~^ for the r^v, /i"*"// and e^v final states, respectively. 

The CLEO-c detector is equipped to measure the momenta and directions of charged 
particles, identify charged hadrons, detect photons and determine with good precision their 



directions and energies. It has been described in more detail previously [7]. Particle identi- 
fication is accomplished using both dE/dx information in the tracking drift chamber and in 
a separate Ring Imaging Cherenkov Detector (RICH) [8]. 

II. DATA SAMPLE AND SIGNAL SELECTION OVERVIEW 

In this study we use 281 pb^^ of CLEO-c data produced in e"'"e^ collisions and recorded 
at the iIj{3770) resonance. At this energy, the event sample consists of a mixture of pure 
D^D^ , D^D , three-fiavor continuum, and '-filj{2S) events. There are also t^t^ pairs, two- 
photon events, and non-DD decays of the ■?/'(3770), whose production rates are small enough 
for them not to contribute background in this study. 

This analysis follows very closely our previous study of D^ -^ fi^u [6, 9]. First we fully 
reconstruct a sample of hadronic D~ decays, that we call tags, and then search for tracks 
that are consistent with a n^ from the decay sequence D^ -^ t^u, r+ -^ n^u, rather than a 
muon directly from two-body D decay. Besides using D^ tags and searching for D^ -^ t^u, 
r"*" -^ Ti^V we also use the charge-conjugate D^ tags and search for D^ -^ '^T^, t^ -^ n^u; 
in the rest of this paper we will not mention the charge-conjugate modes explicitly, but they 
are always used. The loss of rate compared to the muon case, caused by the i3(r+ -^ tt+z/) 
of (11.06±0.11)% [5], is somewhat compensated for by the larger D^ -^ t'^u branching ratio 
as given by Eq. 1. This search has a smeared signal region as compared to the muon case 
because of the extra missing neutrino, and therefore backgrounds are a much more serious 
concern. 

We examine all the recorded events and retain those containing at least one charged D 
candidate in the modes listed in Table I. Track selection, particle identification, 7r°, Ks and 
muon selection cuts are identical to those described in Ref. [9]. 

We have investigated using other r decay modes but they all have significant problems. 
The semileptonic mode euu is embedded in a large D~^ semileptonic background. The p+77 
mode has a MM^ resolution approximately twice as poor, and the vt+tt+tt"!/ mode has several 
additional associated backgrounds, for example D^ -^ tt+tt+tt"??*^ and fifi^u, that severely 
limit its usefulness. 

III. RECONSTRUCTION OF CHARGED D TAGGING MODES 

Tagging modes are fully reconstructed by first evaluating the difference in the energy, AE, 
of the decay products with the beam energy. We require the absolute value of this difference 
to contain 98.8% of the signal events, i.e. to be within approximately 2.5 times the r.m.s 
width of the peak value. The r.m.s. widths vary from 7 MeV in the K^K^n^ mode to 14 
MeV in the K'^n^ii^ii^ mode. For the selected events we then calculate the reconstructed 
D~ beam-constrained mass defined as 



niBC = J^beam " iJ2 fiY^ (3) 

where i runs over all the final state particles. The beam-constrained mass has better resolu- 
tion than merely calculating the invariant mass of the decay products since the beam has a 
small energy spread. 



The m-Bc distributions for all D^ tagging modes considered in this data sample are shown 
in Fig. 2. They are listed in Table I, along with the numbers of signal events and background 
events within the regions shown by the arrows in Fig. 2. The tag candidates are subjected 
to AE and mec cuts explained in our previous paper [6]. The numbers of tagged events 
are determined from fits of the ttibc distributions to a signal function plus a background 
shape. For the background we fit with a shape function analogous to one first used by the 
ARGUS collaboration [10], that has approximately the correct threshold behavior at large 
mBc- To use this function, we first fit it to the data selected by using AE sidebands, defined 
as 5(7 < I Ai^l < 7.5(7, where a is the r.m.s. width of the AE distribution. Doing this mode 
by mode and allowing the normalization to fioat, we fix the shape parameters. For the signal 
we use a lineshape similar to that used for extracting photon signals from electromagnetic 
calorimeters because of the tail towards high mass caused by initial-state radiation [11]. The 
functional form is 
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here A~^ = a. 
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for mec < mo - a ■ am^c 
for rriBc > mo - a ■ arrive 



Here ttibc is the measured mass of each candidate, mo is the "true'' 



(4) 
(or most likely) mass. 



(7mBc i^ ^^^ mass resolution, and n and a are parameters governing the shape of the high 
mass tail. All these quantities are allowed to fioat in the separate fits of each mode. 



Mode 



Signal Background 



K'^TT^TT^ 


77387 ± 281 


1868 


K^TT^TT^TT^ 


24850 ± 214 


12825 


KsTT' 


11162 ±136 


514 


KsTT^TT^TT^ 


18176 ± 255 


8976 


KsTT-TT^ 


20244 ± 170 


5223 


K+K-TT- 


6535± 95 


1271 


Sum 


158354 ± 496 


30677 



TABLE I: Tagging modes and numbers of signal and background events determined from the fits 
shown in Fig. 2. 

We use a total of 158,354±496±475 single tag events for further analysis. The systematic 
error on this number is determined by varying the background function and is estimated at 
0.5%. 



IV. D+ -^ T+u SELECTION CRITERIA 



As in our search for D^ — » fi^u, we calculate the missing mass squared (MM^) defined as 
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FIG. 2: Beam-constrained mass distributions for different fully reconstructed D~ decay candidates 
in the final states: (a) K~^7r~7r~ , (b) K^tt^tt^tt^, (c) Kstt^ , (d) Kstt^tt^tt^, (e) Ks'^~tt^ and (f) 
K^K^TT^ . The solid curves show the sum of signal and background functions. The dashed curves 
indicate the background fits. Events between the arrows are selected for further analysis. 



where Ey^eam is the beam energy, Etrack and {—J^d- — "p^track) ^^^ the measured energy and 
momentum of a single track, assuming that the track is a pion, and 'p'o- is the three- 
momentum of the fully reconstructed D~ . 

The MM^ distribution from Monte Carlo simulations of e+e^ -^ D^D^ , where the D^ is 
fully reconstructed and D~^ -^ t^u, r+ -^ n^V, is shown in Fig. 3. While the pion in this 
decay sequence doesn't have a narrow MM^ peak as in the case of D^ -^ f^^'^, many events 
are in the low MM^ region. The spectrum peaks at low MM^ because the small D^-r^ mass 
difference causes the t'^ to be almost at rest in the laboratory frame and thus the 7r+ has 
relatively large momentum. We must also ensure that we do not accept D~^ — > n~^i> events 
or semileptonic decays with electrons. 
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FIG. 3: Missing mass squared distribution for D^ -^ t^i', t^ -^ tt^u from Monte Carlo simulation. 



Using our D~ event candidates, we search for events with a single additional charged 
track. The crystal calorimeter provides a way of distinguishing this track among muons, 
pions and electrons. We consider three separate cases: (i) the track deposits < 300 MeV in 
the calorimeter, characteristic of a non-interacting pion or a muon; (ii) the track deposits 
> 300 MeV in the calorimeter, characteristic of an interacting pion; (iii) the track satisfies our 
electron selection criteria defined below. Then we separately study the MM^ distributions 
for these three cases. 

We exclude events with more than one additional, opposite-sign charged track in addition 
to the tagged -D, or with extra neutral energy. Specifically, we veto events with extra charged 
tracks arising from the event vertex or having a maximum neutral energy cluster, consistent 
with being a photon, of more than 250 MeV. These cuts are highly effective in reducing 
backgrounds especially from D^ -^ ti^-k^ decays. 

The track candidates are required to be within the barrel region of the detector | cos^| < 



0.81. For cases (i) and (ii) we insist that the track not be identified as a kaon. For electron 
identification we require a match between the momentum measurement in the tracking sys- 
tem and the energy deposited in the Csl calorimeter and the shape of the energy distribution 
among the crystals is consistent with that expected for an electromagnetic shower. 

As demonstrated previously [6], the MM^ distribution has a shape well described by 
two Gaussians for the fi^u mode with a resolution from Monte Carlo simulation (MC) of 
0.0235±0.0004 GeV^. We use different MM^ regions for cases (i) and (ii) defined above. For 
case (i) we define the signal region to be the interval 0.175>MM^ >0.05 GeV^, while for 
case (ii) we define the signal region to be the interval 0.175>MM^ >-0.05 GeV^. Case (i) 
includes 98% of the fi^u signal, so we must exclude the region close to zero MM^, while for 
case (ii) we are specifically selecting pions so the signal region can be larger. The upper limit 
on MM^ is chosen to avoid background from the tail of the K tt"'" peak. The fractions of the 
MM^ range accepted are 46% and 74% for case (i) and (ii), respectively. 

The MM^ distributions for cases (i) and (ii) are shown in Figs. 4 (a) and (b). There are 
12 events in the signal region for case (i) and 8 for case (ii). The electron sample, case (iii), 
shown in Fig. 5, has 3 events in the signal region and is used for background studies. 



V. BACKGROUND EVALUATION 

A. Monte Carlo Estimates 

There are several background sources we need to evaluate. These include background 
from other D'^ modes, background from misidentified D^D events and continuum back- 
ground including that from e^e~ -^ 'jiIj{2S), termed "radiative return." There are a few D~^ 
decay modes that have been identified a priori as possible background sources. These are 
listed in Table II, along with the numerical background estimates we obtain using Monte 
Carlo generation and reconstruction of each specific mode. The branching ratios are from 
the Particle Data Group except for the tt+tt*^ and p+7r° modes where we use new CLEO 
measurements [12]. We note that often at least one photon from the 7r° decay in these two 
modes exceeds our 250 MeV calorimeter energy requirement and causes these decays to be 
vetoed. 

The K 7T^ mode gives a large peak in the MM^ spectrum near 0.25 GeV^. We need to 
evaluate the effects of the tail of the distribution leaking into our signal region. A simulation 
of this background for case (i) and case (ii) yields 2.4 ± 0.5 ± 0.2, and 1.6 ± 0.4 ±0.1 events, 
respectively. The systematic errors are due to uncertainties on the measured branching 
ratios. 

We have also checked the possibility of other D^D~ decay modes producing background 
with an equivalent 1.7 fb~^ Monte Carlo sample; we find no additional events. D^D and 
continuum backgrounds are evaluated by analyzing Monte Carlo samples corresponding to 
4.7 fb~^ and 1.7 fb~^, respectively. To normalize our Monte Carlo events to our data sample, 
we used o'^o-q" = 3-6 ±0.1 nb and fTcontinuum = 14.5 ±1.0 nb [13]. Our total background is 
6.1±0.6±0.3 events in case (i) and 5.0±0.6±0.2 events in case (ii). 
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FIG. 4: The MM^ distributions from data using D^ tags and one additional opposite-sign charged 
track and no extra energetic showers (see text). For the case when the single track (a) deposits 
< 300 MeV of energy in the calorimeter, case (i). The peak near zero is from D^ -^ fi^v events, 
(b) Track deposits > 300 MeV in crystal calorimeter but is not consistent with being an electron, 
case (ii). The arrows indicate the signal regions. The insets show the signal regions with a finer 
binning of 0.002 GeV^. 



B. Background Estimates From Data 



The largest source of background is the tail of the K tt"*" peak. Simulations of the tails 
of distributions, however, are often unreliable. Therefore, we also measure this background 
rate directly from data. 

We select D^D events where one neutral D decays into K^-k^ti^-k^ , K^-k^-k^ or K^-n^. 
These single-tag candidates are reconstructed using tight selection criteria on Ai5 and rriBc- 
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FIG. 5: The MM^ distribution obtained using D^ tags and one additional opposite-sign charged 
track and no extra energetic showers (see text). The track is required to deposit > 300 MeV of 
energy in the calorimeter and be consistent with an electron. The solid curve is data and the dashed 
curve Monte Carlo. 

In this sample, we look for events with only two additional oppositely-signed tracks where 
the RICH system identifies one as a kaon and the other as a pion. We insist that the charge 
of the kaon candidate be opposite to the charge of the kaon in the tag mode. Our aim is to 
isolate the K'^n'^ final state opposite the reconstructed tag signal events. We avoid, however, 
making tight cuts that might ameliorate the effects of tails. 

The K'^n'^ final state is identical kinematically to the K tt"*" state that we wish to emulate 
if we ignore the measurements of the charged kaon and then compute the MM^, as shown in 
Fig. 6 for cases (i) and (ii). 

The event numbers in our signal ranges are 4.8 ± 1.0 ±0.1 for case (i) and 2.5 ± 0.8 ±0.1 
for case (ii). The systematic error arises from the normalization, derived from the fit to 
the MM^ peak near 0.25 GeV^. There are backgrounds, however, in these distributions 
from D -^ 7r"'"7r~ and D -^ ^~tt~^u events where the candidate kaon is a misidentified 
pion. The probability for pions faking kaons in this momentum range has been measured as 
(1.10±0.37)% [8]. Using the known branching ratios for the above two modes, we estimate 
0.08 and 0.17 tt+tt^ events, respectively, and 0.01 and zero 7r+/i~z/ events, respectively, that 
need to be removed from the background estimate, leaving 4.7±1.0 and 2.4±0.7 background 
events. This estimate is in reasonable agreement with the simulation. (Since we are going to 
quote an upper limit in this paper, choosing the Monte Carlo background estimate provides a 
worse limit because less background is subtracted, and thus is the more conservative choice.) 
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Mode 



■S (%) # of events case(i) # of events case(ii) 



vr+vrO 




0.12±0.01 


0.13±0.02±0.01 


1.40±0.07±0.11 






2.77±0.18 


2.44±0.51±0.17 


1.59±0.41±0.11 


yti+Z^ 




0.04±0.01 


1.25±0.03±0.19 


0.46±0.07±0.07 


p+vrO 




0.38±0.03 


0.18±0.05±0.01 


0.23±0.05±0.02 


vrV+i^ 




0.44±0.07 


0.98±0.14±0.15 


0.002±0.001±0.001 


r+i/, r+ ^ 


■ p^v 


0.030±0.005 


0.14±0.01±0.02 


0.15±0.01±0.02 


r+i/, r+ ^ 


■ p^vV 


0.020±0.003 


0.27±0.01±0.04 


0.03(32% C.L.) 


Other D+ 


modes 


- 


0.08(32% C.L.) 


0.08(32% C.L.) 


D^ modes 




- 


0.23±0.12±0.01 


0.42±0.16±0.01 


Continuum 


- 


0.45±0.26±0.03 


0.74±0.33±0.05 


Sum 




- 


6.07±0.60±0.31 


4.99±0.56±0.19 



TABLE IL Monte Carlo estimated backgrounds from all sources. The second errors are systematic 
and are due to uncertainties on the measured branching ratios for D^ background sources and 
production cross-section uncertainties for D^ and continuum sources. The "other D+ modes listed 
at 0.08 at 32% c.l. represent a la upper limit on this contribution. 

Another background check is to both measure the electron background and simulate it. 
We note that the background due to real muons should be almost equal to the background 
due to real electrons. For this study we use the entire MM^ region up to 0.5 GeV^. The 
MM^ distribution due to electron candidates in the data is compared with the one from the 
Monte Carlo in Fig. 5. There are 60±8 electrons in the data compared with 63±3 in the 
Monte Carlo after normalizing to the luminosity in the data. (In the signal region there are 
3 events in the data versus 3.9±0.1 in the Monte Carlo.) The good agreement establishes 
that the Monte Carlo properly predicts the semileptonic decay backgrounds. 

VI. BRANCHING RATIO LIMITS 



We do not observe a statistically significant difference between the number of signal and 
background events. For case (i) we have a net signal of 5.9 t^v, r+ — * -k^V signal events and 
for case (ii) our yield is 3.0 events For each of our two cases denoted by j, where j represents 
either case (i) or case (ii), the expected number of events 
Brv = B{D^ -^ T^v) through the relationship 



^expected' is related to the true 



expected 



-'^tags X ^T 



^B{t^ 



TT 



V) X e^ + nI 



bkg' 



(6) 



where A^tags is the number of single tag events and equals 160,729, after correcting for the 
slight difference in reconstruction efficiency for tags opposite a single track versus tags op- 
posite a typical D^ decay; e^ is the efficiency, and N^^^ is the background. The Poisson 
probability distribution Lj{Bru) for B{D^ -^ t^v) in each case j has a mean equal to 



N\ 



expected 



and is given by: 



Lj{Br,) = (^^j X exp (-iV^p^^^^d) X (^e'xpected)^ 



(7) 
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FIG. 6: The MM^ distribution from data events with a single D^ or D tag and the other neutral 
D decaying into two tracks, most likely K^tt , where the kaon information is ignored. For the two 
cases: (a) track deposits < 300 MeV of energy in the crystal calorimeter and (b) track deposits 
> 300 MeV in the calorimeter. The arrows delineate the relevant signal regions. 



where A^-' is the number of detected r^u candidates: 8 for case (i) and 12 for case (ii). 

Our results for the branching fractions are found by doing a simultaneous likelihood fit 
of the distributions described in Eq. 7. We take into account the the different efficiencies in 
cases (i) and (ii) that arises from both the MM^ acceptance (46% and 74%) and the efficiency 
of not having another unmatched shower in the event with energy greater than 250 MeV 
(93.9% and 91.8%). We have previously found [6] that the Monte Carlo matched within 
1.8% our measurement of the extra unmatched shower cut and thus use a slightly larger 2% 
for the systematic error on this quantity. Overall, the efficiencies are 18.7% (22.4%), for case 
(i) and (case(ii)), respectively. 



We find B{D^ 



T^V) 



[l.Stll ± 0.1) X 10-3 and (0.8+^1 ± 0.2) x 10-^ for cases (i 



and (ii), respectively, where the statistical errors result from the values of Bru corresponding 
to 34% of the area under the Lj distribution above and below the maximum value. 

The errors on the backgrounds are treated as systematic and are obtained by varying 
the background contributions in the likelihood distribution. The systematic errors on the 
branching ratio from sources other than backgrounds are listed in Table III; they are neg- 
ligible in comparison with the statistical uncertainty. A more detailed explanation of the 
sources of systematic errors can be found in our previous Letter [6]. 
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Systematic errors (%) 



MC statistics 


0.2 


Track finding 


0.7 


PID cut 


1.0 


Minimum ionization cut 


1.0 


Number of tags 


0.5 


Extra showers cut 


2.0 


Total 


2.6 



TABLE III: Systematic errors on the D^ -^ t'^u-t branching ratio. 

To obtain a combined result for the branching fraction, we construct the global likelihood 
as the product of the two Poisson probability distributions, and we extract the value of Erv 
which maximizes this likelihood function. We find B{D^ — > t^v) = {1.2'to'l±0.l) x 10~^. We 
caution the reader that this is not a definitive measurement but an intermediate step used in 
the process of forming an upper limit. (Had we used the data to estimate the background, 
the branching fraction would be lower.) 

Since the result is not statistically significant we quote an upper limit of 

i3(D+ ^r+z/) < 2.1 X 10-^ (8) 

at 90% confidence level. 

The ratio to the expected rate in the Standard Model using our measured B{D^ -^ f^^^) 
is < 1.8 at 90% confidence level. 



VII. CONCLUSIONS 

We have measured the first upper limit on the decay D^ -^ r^u. We limit B{D^ — > r^u) 
branching ratio to < 2.1 x 10~^ at 90% confidence level. We use our previously measured 
result of B{D+ -^ /i+z/^) = (4.40 ± 0.66^^^) x 10-^ [6], coupled with the evaluation of 
Eq. (2) of 2.65, to limit the ratio R = T {D+ -^ t+v) /V {D+ -^ fi+u) to that expected in 
the Standard Model to < 1.8 at 90% confidence level. Thus lepton universality in purely 
leptonic D~^ decays is satisfied at the level of current experimental accuracy. 
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